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We discuss a specific entanglement distillation scheme under the constraint of finite samples of 
entangled qubit pairs. It is shown that an iterative process can be explicitly formulated. The average 
fidelity of this process can be enhanced by introducing conditional storing of entangled qubit pairs 
in each step of the iteration. We investigate the corresponding limitations on the size and the initial 
fidelity of the sample. 
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I. INTRODUCTION 

The quantum concept of entanglement is the most in- 
triguing feature that allows to establish new physical 
paradigms in information processing. The corresponding 
non-classical protocols need a certain degree of entangle- 
ment as a quantum resource 0]. Consequently one has 
to understand how entanglement can be processed and 
how it can be measured (2 HI • 

In this respect the typical situation is the following. 
Several parties share components of an entangled sys- 
tem. Processing of entanglement then means that they 
can perform local unitary operations and local measure- 
ments on their respective parts of the complete system. 
Under these conditions it has been shown 0. S IE 0- 
that two parties, supplied with non-maximally entangled 
qubit pairs, can extract a sample of stronger entangled 
pairs. But so far no general approach to optimal distil- 
lation exists. However, in order to improve the fidelity, 
modifications of the original versions have been proposed 
0, El 0|. Lower bounds for the fidelity of entangle- 
ment distillation based on faulty local operations have 
also been studied (l^ . 

In the present work we investigate the entanglement 
distillation protocol described in [6] for the case that only 
a finite sample of entangled qubit pairs is available. In 
particular, we propose a relatively simple, iterative distil- 
lation scheme that starts from a finite number of identical 
pairs and delivers a distilled pair applicable for further 
communication tasks. The behaviour of the correspond- 
ing mean fidelity turns out to be particularly interesting 
for small initial samples of entangled qubit pairs. 



II. ENTANGLEMENT DISTILLATION 

For any non-classical communication, entangled sys- 
tems first have to be distributed between a sender (Al- 
ice) and a receiver (Bob). In the course of this distribu- 
tion the systems are influenced by various noise sources, 
which reduce the amount of entanglement. Any distilla- 



tion process for such distributed systems is restricted to 
local operations and classical communication (LOCC) of 
the parties. Moreover, realistic entanglement distillation 
schemes have to take into account errors 0] and the fact 
that Alice and Bob share only a finite amount N of mixed 
entangled systems. In the present paper we focus on the 
latter restriction regarding resources. 

The explicit distillation protocol we refer to was intro- 
duced in Refs. 0,@. In the present work we only study 
this specific distillation process under the constraint of fi- 
nite resources. However, other highly effective distillation 
protocols have been proposed in the literature. Their ex- 
clusion in the present discussion does not mean that they 
cannot be applied to finite resources. In particular, the 
quantum hashing protocol @ can be well applicable to 
finite resources and presumably leads to highly distilled 
qubit pairs. It will, however, need a considerable effort 
for the corresponding book keeping of the so-called likely 
sets. 

The process jj, |g] conditionally increases the fidelity 



F(p) = tr (* + /5) = A 



of the Bell-diagonal mixed state 



A<f> 



(1) 



(2) 



* E-mail me at: Stefan.Probst@uni-ulm.de 



describing a qubit pair, where we introduced the abbre- 
viations * ± = \4> ± )(4> ± \ and = \ip ± ){4> ± \ for the four 
Bell-states. 

This state is non-separable if any of the coefficients 
A, B, C or D is larger than 1/2 Q. Without loss of 
generality we choose A > 1/2. 

The central element of the protocol is a CNOT trans- 
formation. We briefly recall the basic ideas 0,0]. Two 
qubit pairs 1 and 2 each described by the state p, Eq. (2J, 
are processed in one step. Alice holds the qubits 1a and 
2a and Bob holds the qubits Is and 2b- These qubits can 
now be treated locally using a sequence of operations 
(I) Alice and Bob rotate their qubits locally. These rota- 
tions exchange the \&~ contribution with the <& _ contri- 
bution of the initial state p, Eq. Q [2fl| - (II) Alice and 
Bob then perform CNOT operations on their respective 
qubits. The qubits of pair 1 (qubits 1a and 1b) act as 
control qubits. (Ill) Both measure their target qubits 
(qubits 2a and 2b) in the computational basis {|0), |1)}. 
They obtain either the result »0« or »1«. (IV) Alice and 
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Bob classically communicate their measurement results. 
The distillation is successful with probability 

p^ = (A + B) 2 + (C + D) 2 (3) 

if the combined result reads »00« or Then they 

keep pair 1 now described by the conditioned density 
operator 



1 



(A 2 +B 2 )®+ +2CD*" 



+ (C 2 + D 2 ) *+ + 2AB * 



(4) 



If Alice and Bob read off the results »01« or »10«, the 
reduced density operator of pair 1 reads 

1 



1 - p( s ) 

f (AC + BD)^ 



(AC + BD)<& + + (AD + BC)V 



(AD + BC)^- 



(5) 



In this unsuccessful case Alice and Bob discard pair 1. 

In the successful case pair 1 is mapped from a Bell- 
diagonal state p, Eq. to the Bell-diagonal state p( s \ 
Eq. The corresponding new fidelity is given by 



F(p^) = 



A 2 + B 2 



(A + Bf + (C + Df 



(6) 



which is higher than the initial fidelity for A > 0.5. 

However, this fidelity does not completely describe the 
single-step process. A complete description also has to 
take into account the unsuccessful case. The distillation 
process then leads to the conditioned density operator 
p( u \ Eq. <HjJ|. Hence if we are interested in the perfor- 
mance of this particular distillation process, we have to 
take into account the fidelity F(pW) < \. On the other 
hand, using local operations and classical communication 
the parties can always generate two qubits with fidelity \ 
if the unsuccessful case occurs. Applying such additional 
operations we can define the average fidelity 



(F)=p^F(p^) + (l-p^ 
= A + B(1-2A) 



(7) 



for a single step process. For all possible values of B the 
average fidelity (F) turns out to be equal to or less than 
the original fidelity A. This is of course also true for the 
average fidelity 

(P)=pUF(pM) + (l-pM)F^) 
= A 2 + A(C-B) + B(1-C) 

based on the density operator p( u \ Eq. ©. This is con- 
sistent with the fact that the total amount of entangle- 
ment cannot increase under a LOCC process 0]. 

However, we emphasise that so far we have just dis- 
cussed a single distillation step performed on two entan- 
gled qubit pairs. The situation becomes different and 
more interesting when we now consider a finite ensemble 
with N > 2 pairs. 



III. DISTILLATION FOR A FINITE SET OF 
ENTANGLED SYSTEMS 

Alice and Bob now perform the distillation scheme 
with an even number TV of pairs [2j] ■ 

After one distillation step they have j < N/2 pairs left 
with probability 



p(N,j) = 



(T 



p 



(s)]N/2-j 



(9) 



which contains the success probability p( s \ Eq. (0). 

The totally unsuccessful case occurs when all pairs 
have to be discarded, that is j = 0, in a single distil- 
lation step. Hence starting from N pairs we can define 
the average fidelity 

(F)(N) = P (N,0) F^ u) +[l-p(N,0)]F(p^) (10) 

using the density operators p^ s \ Eq. Q, and the fi- 
delity F( u ) for the unsuccessful case. With an increasing 
amount N of pairs the average fidelity (F)(N) also in- 
creases because of the decreasing probability p(N, 0) to 
lose all qubit pairs. Using Eq. (fflTTfl one finds that 
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N min = 2- 



A — F(p^) 

JP(n) -F(pW) 



ln(l — p 



(11) 



pairs are needed to obtain an average fidelity 
(F)(Af m j n ) = A. This minimal number A m ; n of qubit 
pairs shows a strong difference depending on the choice 
for ^( u ). Fi gure shows the minimal size Amin for a 
finite ensemble described by Werner-states (B — C — 
D = (l-A)/Z) 0. If we simply choose F^ = F(p^) 
to characterise the CNOT distillation itself we obtain a 
strong dependence (dashed curve) on the initial fidelity 
A. In particular, A^in diverges like ln(vl— \) for A — * 
If, on the other hand, we substitute the LOCC bound- 
ary 

F (u) = 1 

, we always get finite values for Amin which 
only depend weakly on A. In the plot one can therefore 
clearly identify the region where on average the fidelity 
increases, that is (F)(N) > A. 

Moreover, if j > 2 pairs are left after such a first step 
we can continue with the distillation. We consider such 
an iteration in the following paragraph. 



IV. ITERATIVE DISTILLATION 

The resulting state p^ of a successful distillation step 
is again Bell-diagonal and hence the complete process can 
be applied iteratively, as long as qubit pairs are left to 
use. Starting from an initial density operator p^ = p, 
Eq. J2J, the density operator pf} 1 is mapped on a density 

operator fff 1 if the i th distillation step was successful. 
The corresponding coefficients of the Bell-projectors, see 
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FIG. 1: Minimal size iV m j n , Eq. jlip . of the finite sample of 
entangled qubits to obtain an average fidelity (F) (JV m ; n ) = A 
in a single step distillation. The curves have been calcu- 
lated for specific Bell-diagonal states (Werner-states) with 
B — C — D — (1 — j4)/3. Nevertheless, they show the generic 
behaviour. The dashed curve results from the substitution 
p(u) _ p(pw] m Eq. ITU and diverges for small initial fi- 
delities. In contrast to this, the minimal size for the LOCC 
choice = | stays always finite (solid curve). Clearly for 
4^1 we obtain iV m i n = 2. 



Eq. J3J, transform as 
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(12) 



This 



with the corresponding success probability 
mapping was studied in detail in Ref. [l^j . 

If, however, the i th step was unsuccessful we have in 
principle a mapping from pf} 1 to In general, the 

density operator , see also Eq. (§J , has a lower fidelity 
than g , which is the fidelity of the density operator that 
can be generated by a LOCC process. If the step was 
unsuccessful we therefore assume from now on that Alice 
and Bob perform appropriate operations to prepare local 
qubits with fidelity = \. 

A full iteration of the entanglement distillation for fi- 
nite quantum resources has to take into account all pos- 
sible combinations of these successful and unsuccessful 
steps. Our aim then is to obtain at the end of a CNOT 
distillation an entangled qubit pair with an average fi- 
delity as high as possible. 

Note that one can also ask a different question for a 
small finite sample. Can we distil a single qubit pair with 
very high fidelity if we allow for a lower success proba- 
bility? That is, we do not care if the distillation fails 
several times, but when it is successful it must deliver a 
highly entangled pair. Then the CNOT distillation pro- 
cess seems to be not very suitable. In this case quantum 




Save one pair 
n := n — \ 



Purify n qubit pairs 
j purified pairs with 
probability p(n,j),Eq. (9) 




Successful distillation 
j > : Fit = Fidelity of pairs of last iteration 
j = : fit = Fidelity of last saved pair 




Unsuccessful distillation 

fit ~ \ 



FIG. 2: Flow chart for the iterative CNOT distillation scheme 
with finite quantum resources. Starting from N qubit pairs 
we extract a distilled pair. In every step n/2 projective mea- 
surements are performed and some qubit pairs have to be dis- 
carded. To increase the average fidelity a backup pair can be 
saved in every iteration step for an odd number n. Moreover, 
the iteration stops if only two pairs are left after the projec- 
tive measurement and no backup pair was saved previously. 
Rarely the distillation may still fail if no pair was saved and 
all pairs after the last projection have to be dismissed. The 
scheme also defines the corresponding fidelity Fa of the iter- 
ative process. Averaging over many runs we obtain (Fi t )(N) 
for N qubit pairs. 



hashing can turn out to be a powerful method. 

The complete iterative CNOT distillation process for 
finite resources can be expressed in algorithmic form. For 
a given initial density operator the process depends only 
on the initial amount of entangled qubit pairs shared by 
Alice and Bob. An important iterative improvement of 
the process can be achieved for odd numbers of qubit 
pairs. In each step of the iteration we may obtain an odd 
number n of pairs and hence in the simplest case we can 
store the additional pair as a backup [ifj. This backup 
pair can be used whenever we would have to discard all 
pairs in some further step of the iteration. 

The stop condition of the algorithm depends on the 
specific aim of the process. As we have already seen in 
Eq. , it is impossible to obtain an increase of the aver- 
age fidelity for two qubit pairs. So if we end up with only 
two pairs left after a distillation, it seems to be unwise 
to continue with the iteration, in particular if we have 
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stored no backup so far. The corresponding flow chart 
of such an algorithm is shown in Fig. [21 The average fi- 
delity (F it )(N) of the iterative process can be computed 
using this algorithm. 

In the following example we present the behaviour of 
such an iterative distillation process. We simulate the 
process using the probability p(n,j), Eq. ©, to obtain 
j qubit pairs when starting from n. Moreover, we need 
the mapping of Eq. l(T2ll . 



V. EXAMPLE OF THE ITERATIVE 
DISTILLATION FOR SMALL FINITE SETS 
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To demonstrate the behaviour of the iterative distilla- 
tion we have again chosen Werner-states (B^ = Cq^ = 

= (1 - A ( s) )/3) as initial states. 

(Fit)(N)/4 s) 



1.10 -- 




FIG. 3: Relative average fidelity (F it )(N) /A ( g' of the iteration 
as a function of the initial fidelity A^ for different amounts 
N of qubit pairs. The initial states are again Werner-states 
(B ( S) = C ( S) = D ( S) = (1 - 4 S) )/3). For N = 4 qubit pairs 
(solid line) we only get an improvement if Aq > 0.65. For 
greater amounts of qubit pairs we always get an improvement. 
The advantage of the backup can be seen by comparing the 
two cases N — 5 (dotted line) and N = 6 (dashed line). Due 
to the fact that for odd numbers we always have a backup 
pair, the average fidelity becomes higher as compared to the 
surrounding even cases. 

First we compare in Fig.Elthe relative average fidelity 
(F lt )(N)/A^ as a function of the initial fidelity for 
small resources of qubit pairs. In accordance with the 
algorithm of the previous paragraph we have calculated 
the final average fidelity (F lt )(N) with backup for N = A 
(solid line), N — 5 (dotted line) and N = 6 (dashed line) 
pairs. We see that N = 5 pairs are on average always 
superior to the nearby even cases, because for an odd N 
already in the first step of the distillation a backup pair 
is stored. For N — 4 there is a minimal initial fidelity 
which is needed to have a successful iterative distillation. 

Second, we show in Fig. the generic behaviour of 
the average fidelity (F it )(N) depending on the avail- 
able amount N of pairs. We compare the iterative 



FIG. 4: Comparison of the average fidelities of an iteration 
depending on the amount N of initially accessible qubit pairs. 
The plot shows the average fidelities of iterative schemes 
without backup (dotted curve), with backup (dashed curve), 
FigEl an d the mapping of the completely successful case (solid 
curve), Eq. Ijl2fl. We show their behaviour for the mentioned 
example state with an initial fidelity Aq — 0.75. The plot 
again shows the advantage of storing one pair in the case of 
odd numbers JV. It turns out that for even numbers N it is 
on average advantageous to drop one pair and to start the 
distillation with one pair less. This shows that simple backup 
schemes improve the CNOT distillation and more sophisti- 
cated methods can certainly be discussed. 



results without and with backup for a fixed initial fi- 
delity ^4q S ^ = 0.75 to the completely successful mapping, 
Eq. 11211 . Storing a backup pair in the case of an odd 
amount clearly leads to a higher average fidelity. The 
zigzag in the curves reveals the difference between odd 
and even numbers N of pairs. If the iteration starts with 
an odd number of entangled qubit pairs, at least one pair 
is always saved. One can even see that for even numbers 
it is on average better to drop one pair before performing 
the first distillation step. Although this means to start 
with one pair less, the average fidelity is higher. This is 
most obvious for small finite ensembles. 



VI. CONCLUSIONS 

In the present paper we have analysed an iterative 
scheme for a known distillation protocol. This protocol 
is especially useful for the implementation of an iteration 
since it puts very low restrictions on the density oper- 
ators of the processed systems. In particular, we have 
emphasised the application to finite quantum resources. 
In this respect it was possible to demonstrate the limi- 
tations on the needed number of entangled qubit pairs 
as well as on their initial entanglement for a successful 
distillation. Stronger entanglement can be obtained it- 
eratively already for small initial numbers of pairs, even 
though many pairs have to be sacrificed in order to obtain 
one distilled pair. We have lowered this loss by introduc- 



5 



ing backup pairs in our algorithmic description. 

Clearly this is not the only possibility to achieve itera- 
tive entanglement distillation. First, one can think of re- 
cycling backup systems in the distillation process. This, 
however, works only in a rather narrow regime in which 
the fidelities of the recycled pairs are close. Second, it 
is also possible to iterate completely different distillation 
methods, like quantum hashing |4( for finite resources. 
Both possibilities need a considerable amount of book 
keeping in contrast to the iteration presented here. 

Finally, it would be important to simulate distillation 
for an experimental system, which offers a way to control 



a finite number of qubit pairs. A possible candidate for 
this would be an optical lattice filled with atoms, which 
can be controlled by collisions [Til ITii. Hfll] . 
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